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RESUME 


On  a  concu  un  module  math^matique  de  reticule  de  t&te  chercheuse 
h  modulation  d’amplitude  et  d^termin^  sa  position  d*6auilibre  stable  en 
presence  d*une  cible  b  double  points,  ou  h  simple  point  accompagn^e 
d'une  cible  ^tendue.  Le  comportement  de  ce  reticule  constitue  I’objet 
d'une  simulation  sur  ordinateur  num^rique.  (NC) 


ABSTRACT 

A  generic  AM-reticle  seeker  was  modelled,  and  its  steady-state 
equilibrium  position,  in  the  presence  of  either  two  point  tarqets  or  a 
point  and  an  extended  target,  determined  by  digital  computer 
simulation.  (U) 
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1 . 0  INTRODUCTION 

An  infrared  seeking  missile  employs  the  radiation  from  the  target 
as  its  homing  source;  its  qimballed  seeker  head  nrecesses  to  follow  this 
source.  Although  there  has  been  a  rapid  development  of  charge-couoled 
devices  suitable  for  processing  target  radiation,  current  threats  are 
more  likely  to  use  mechanical  scanning.  The  mechanical  scanners,  which 
have  fields  of  view  of  a  little  more  than  a  degree  to  several  degrees, 
are  either  of  the  rosette  or  reticle  types.  The  rosette-type  scans  with 
a  very  narrow  instantaneous  field  of  view,  thus  permitting  the 
recognition  of  individual  characteristics.  The  reticle  scanners  include 
those  which  provide  pulse-position  modulation,  freouency  modulation  or 
amplitude  modulation  to  encode  the  angular  position  of  the  target 
relative  to  the  seeker  axis.  In  this  document  a  seeker  employing  an 
amplitude-modulating  (AM)  reticle  is  examined;  related  reticles  are 
employed  in  certain  American  and  Soviet  missiles. 

With  two  targets  within  its  field  of  view,  the  seeker  of  a 
missile  employing  an  AM  reticle  has  a  tendency  to  point  nearer  to  the 
more  intense  one,  a  property  which  is  employed  to  decrease  the 
effectiveness  of  the  infrared  seeking  missile;  a  taroet  can  dispense  a 
decoy  flare,  which,  if  it  is  sufficiently  intense,  may  be  accented  bv 
the  missile  seeker  as  a  preferential  source  of  radiation  to  the  target 
itself.  Analyses  of  infrared  decoy  flare  effectiveness  require  the 
determination  of  the  position  to  which  the  seeker  of  the  at tackina 
missile  will  point  relative  to  the  tarnet  and  flare  positions.  Tn 
simplify  the  calculations,  it  is  often  assumed  that  the  seeker  will 
point  to  the  weighted  mean,  that  is:  the  product  of  image  intensity  and 
distance  from  the  centre  of  the  reticle  will  in  the  steady  state  be 
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equal  for  the  target  and  the  flare.  This  assumption  may  be  valid  when 
two  infrared  sources  are  far  from  the  missile,  and  the  images  on  the 
reticle  are  thus  small.  But  it  will  no  longer  be  true  as  the  missile 
approaches  the  targets,  and  the  images  occupy  a  substantial  portion  of 
the  reticle. 

A  model  of  a  generic  AM  seeker  is  developed,  and  its  performance 
examined  for  certain  specific  cases  which  are  tractable  by  exact 
mathematical  solutions.  As  well  as  introducing  the  reader  to  the  seeker 
behaviour,  these  examples  serve  to  validate  the  computer  programme.  A 
digital  programme  is  developed  from  the  model,  and  employed  to  determine 
the  seeker  equilibrium  in  the  presence  of  two  point  taroets  as  a 
function  of  target  separation  and  relative  target  intensities.  This  is 
then  extended  to  include  the  case  where  one  target  is  an  extended 
source.  The  results  nay  be  of  assistance  in  the  modelling  reouired  for 
flare-effectiveness  studies. 

This  report  covers  the  work  performed  at  DREV  between  1  October 
1979  and  1  August  1980  under  PCN  21313,  '‘Multiple  Target  Tracking  with  a 
Reticle  Seeker”. 
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2-0  MODEL 

2.1  General 

A  generalized  schematic  of  an  infrared  seeker  is  shown  in  Fin. 

The  infrared  emission  from  the  target  is  focussed  onto  a  reticle 

sninning  about  its  centre  at  a  freouency  of  f^.  The  reticle  consists  of 

alternate  opaque  and  transparent  sectors,  which  chop  the  incident 

radiation  to  produce  a  carrier  of  frequency  p.f  at  the  photocell 

s 

output,  where  p  is  the  number  of  pairs  of  opaque  and  transparent 
sectors.  The  chopping  efficiency  is  maximum  for  an  image  capable  of 
being  instantaneously  masked  or  unmasked  totally;  such  image  can  be 
either  a  point  or  a  straight  line  (or  its  extension)  passino  throuqh  the 
centre  of  the  reticle.  This  waveform  is  modulated  at  the  spin  frequency 
by  a  semi-transparent  phasing  sector.  An  automatic  qain  control  (AGC) 
keeps  the  amplitude  of  the  carrier  constant  and  independent  of  the 
target  intensity,  which  is  a  requirement  for  linear  control.  This  signal 
is  demodulated  and  produces  quadrature  components  at  the  spin  freouencv, 
the  amplitudes  of  which  are  a  measure  of  the  angular  error  between  the 
look  angle  of  the  seeker  and  a  line  to  the  target.  These  error  signals 
are  applied,  via  torauinq  motors,  to  nrecess  the  reticle  assembly  to 
reduce  the  error  angle  to  a  small  value. 

2 . 2  Reticle 

The  reticle  used  in  this  model  is  shown  in  Fin.  2.  It  possesses 
twelve  pairs  of  transparent  and  opague  spokes,  and  the  inner  boundaries 
of  the  phasing  sector  are  delimited  by  semi -circles ,  centred  on  opnosite 
sides  of  a  diameter  of  the  reticle,  the  radiuses  of  which  are  half  that 
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FIGURE  1  -  Control  system 


The  reticle 
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of  the  reticle.  The  phasing  sector  has  a  transmission  of  so  that 

the  average  signal  will  be  the  same  during  image  dwell  nn  the  phasinn 
sector  and  image  dwell  on  the  spokes;  this  will  tend  to  improve 
discrimination  against  large  images.  Harmonic  reference  signals  are 
generated  at  the  spin  freouency  of  the  reticle,  which  rotates  in  a 
clockwise  direction  at  a  rate  of  rad /s . 

To  illustrate  the  manner  in  which  the  reticle  encodinn  is  used  to 
produce  error  signals,  an  example  is  now  given  in  which,  for  ease  of 
comprehension,  the  signal  processing  is  greatly  simolified.  Let  the 
target  image  be  at  a  distance  r  from  the  centre  of  the  reticle.  Then, 
from  the  reticle  geometry  given  in  Fig.  3,  it  can  be  seen  that  the  image 
will  be  interrupted  by  the  phasing  sector  for  an  angle  of  2  arcsin  r(=2R) 
for  each  rotation  of  the  reticle.  For  this  narticular  figure  only,  the 
image  has  been  rotated  around  the  centre  of  the  stationary  reticle. 


FIGURE  3  -  Reticle  geometry 
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Suppose  that  the  image  is  at  shown  in  Fin.  ^(a),  which 
has  been  drawn  for  the  instant  that  the  reticle  and  spatial  references 
are  coincident.  The  photocell  sional  will  then  he  as  in  Fig.  ^(h) 
Reference  signals  (Fig.  4(c))  are  generated  by  the  reticle  motor,  and 
are  phased  so  that  the  information  contained  in  the  photocell  signal  can 
be  correctly  recovered.  The  photocell  signal  is  rectified  (Fig.  4(d); 
and,  in  conjunction  with  the  reference  signals,  its  5INF  and  COSINE 
components  are  determined. 


COSINE  = 


1  r  271- 3-4;  2tt 

—  r cos  a  da  +  f  cos  a 

""  Lo 


—  sin  B  cos  \p  , 

7T 


[1] 


where 


=  -  —  r  cos  , 


a  =  u)  t  . 
s 


Similarly 


SINE  =  -  r  sin  4^ 


[2] 


The  spatial  axes  are  defined  such  that  when  4;  z  0,  the  imane  lies  on  the 
positive  X  axis,  and  when  \i>  r  tt/2,  it  lies  on  the  positive  Y  axis.  It 
is  required  that,  when  an  error  is  in  the  positive  direction,  the  error 
signal  be  also  positive.  For  this  to  occur,  the  error  signals  for  the 
X  and  Y  directions  (E^  and  E^)  must  be  equal  to  (-COSINE)  and  SINE 
respectively,  which  requires  the  use  of  ( -cosine )  and  sine  reference 
signals.  Flenceforth,  these  will  be  the  reference  sinnals  employed. 
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TARGET 

IMAGE 
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(  d  ) 


FIGURE  4  -  Waveforms  required  for  example 
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As  this  simplified  example  demonstrates  a  linear  control  it  would 
be  anticipated  that  the  more  sophisticated  model,  which  is  emoloyed 
throughout  the  remaining  analyses,  will  show  a  similar  tendency, 

2.3  Signal  Processing 

The  photocell  signal  is  first  filtered.  In  the  model,  only  the 
carrier  (at  the  spoke  frequency)  and  its  immediate  side-bands  are  oassed 
by  the  filter.  If  the  spin  frequency  is  f  ,  that  of  the  carrier  will  be 

o 

12f  (assuming  a  12— spoke  reticle),  and  the  filter  will  pass  only 
s 

components  at  i2f  and  13f  .  In  practice,  the  signals  from  the 

s  s  s 

various  target  images  are  summed  by  the  photocell,  and  then  filtered. 
As  the  signal  processing  to  this  point  is  linear,  the  model  achieves  the 
same  result  by  filtering  the  signals  from  the  individual  images,  and 
then  summing.  Thus  the  filtered  signal  from  the  i*th  inane  ran  be 
represented  by 

13 

'’ij  •  [3j 

The  total  filtered  signal  is  then 

S(F)  =  I  S(F  )  .  [4] 

X 

1 

Before  this  signal  is  detected,  some  form  of  automatic  nain 
control  must  be  applied,  otherwise  the  amplitude  of  the  error  sinnals 
will  be  a  function  of  the  image  intensities,  and  linear  control  would 
not  be  feasible.  The  type  of  AGC  employed  herein,  is  that  of  dividing 
5(F)  by  its  root-mean-square  (RMS)  value,  to  give 


S(AGC)  =  S(F)/RMS, 


[5] 
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where 


13  3  3  i 

RMS  =  [  Z  (a;=  +  b^)  ] 

j=ll  ^  ^ 


[6] 


and 


a.=Za..;  b.  =  2^b..  . 


j  i 


j  i  11 


This  sinnal  is  then  square-law  detected,  and  the  coefficients  of  sin 


and  cos  a  extracted.  Let 


13 


S(AGC)  =  I  (c .  sin  ja  +  d,  cos  ja)  . 


j=ll 


3 


3 


[7] 


Then 


S'^(AGC)  =  cos  a  +  Ey  sin  a  +  higher  frequency  components,  [8] 


where 



Ey  =  ■  ^15^  ‘ 


The  hiqher  frequency  conoonents  are  filtered  out,  and  and  fy  ^re  used 
to  drive  the  reticle  assembly,  in  the  X  and  Y  directions  resnect ivelv , 
via  the  toroue  motors  which  have  a  gain  of  G. 


A  flowchart  of  the  model  is  given  in  Fiq.  5. 

2.4  Error  Signal  for  Small  Error 

If  the  target  image  is  displaced  only  a  small  distance  from  the 
reticle  centre,  then  sine  and  cosine  approximations  can  be  emnloyed  to 
simplify  the  calculations.  This  will  permit  a  further  examnle  of  the 
signal  orocessinq  to  be  given,  without  as  yet  resortino  to  a  computer. 
In  Fig.  6(a)  is  shown  the  photocell  signal  for  a  tarnct  imane  on  the  X 


A 
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axis  at  a  distance  r  from  the  reticle  centre^  where  r  is  nuch  s^^aller 
than  ff/p,  and  p  is  the  number  of  pairs  of  sectors  (=12). 


6  e 

!!♦  Hh* 

(d) 


nnniiminji^^ 


FIGURE  6  -  The  photocell  waveform 


The  filtered  signal  then  has  the  following  form 

1 

S(F)  =  Z  [a  ^^sin(p  +  i)a  +  ^  t9] 

i=-l  ^  ^ 


The  equation  will  now  be  written  in  the  following  tabular  form 


S(F)  = 


^“1  ^1 


-1 


V  V 
0  1 


[10] 


where  u  is  the  coefficient  of  sin(n  +  p)o, 
n 

and  is  the  coefficient  of  cos(n  +  p)a;  n  =  -1,  0,  1, 

and  the  straight -line  brackets  are  a  f unction  which  performs  the 

following  operation 




i 


KNCLASSIFKD 

1> 


1 

Z  sin(n  +  p)  a  +  cos  (n  +  p3a. 


The  waveform  of  Fig.  6(a)  will  be  analysed  by  first  separatinq  it  into 
the  two  parts  shown  in  Fiq.  6(b).  It  follows  that  the  5(F)  for  the 
photocell  signal  is 


S(F) 


1  -  cosCp  “  1)^ 

(p  -  l)7r 
0 


-2  1  -  cos  p3 

—  + - 

7T  pTT 

0 


1  -  cosCp"**l)B 

(P  +  l)ir 
0 


[11] 


Assume  nr<<l,  so  that  6  =  (arcein  r)  ==  r,  then  eo.  11  can  he  written 


S(F)  = 


(P  -  1) 

Zrr 

0 


~  2  (p  +  1) 

IT  2v 

0  0 


[12] 


Let  the  target  image  be  displaced  from  C(x,0),  or  P(r,0),  to  C(x,y),  or 
^(^,4^),  where 

r  =  (x  +  y  )*  , 

\\f  =  Arctan  y/x  . 


In  the  model  the  phase  angle  is  introduced  after  the  filtering,  so  that 


2 

(p-lj  r  cosCp-1)^ 

2  TT 


cos  p4^ 

7r 


2 

(p-*-!)  r  cos(p+l)i() 
2  TT 


S(F)  = 


(P-l)  r  sin(p-l)i|/ 
2  TT 


sin  pi;/ 

tr 


(PID  sin(p+l)(// 
2  TT 


In  the  case  of  multiple  targets,  it  is  at  this  point  that  the  S(F.) 
would  be  summed  to  form  S(F).  The  RMS  value,  because  r  is  small,  is 
very  close  to  2/tt  .  Hence 


- 
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(p-I)  r  cos(jp“l)i^ 
4 

2 

(P-I)  r  sin(p-l)i(/ 
4 


-cos  plj/ 
-sin  pi|> 


2 

(£M)  r  cosCp+1)^ 
4 

(P’^’l)  r^sin(p+l)i|; 
4 


[14] 


Equation  14  is  substituted  in  ea.  8  to  give 


{cos  pi|>[  (p-l)cos(p-l);|;  +  (p+3)cosCp+l)\|^] 

4 

+  sin  pip  [  (p'-l)sin(p-l)i^  +  (p+l)sin(p+l)\(»] } 
2 

=  6  r  cosi^ 


2  r  r 

c  -  {cos  p4/  [  (p-l)sin(p-l)i|;  -  (p+l)sin(p+l)\(>] 

Y  4 

-  sin  [  (p-l)cos(p-l)*  -  Cp+l)cos(p+l)i/;]} 

-  ^  2  .  . 

=  6  r  sinij^ 


[15] 


[16] 


It  is  seen  that  for  small  disolacements  the  system  has  a  snuare-lpw 
control . 

2.5  Equilibrium  Position  for  Two  Targets  (Small  Displacement) 

A  further  example  is  now  given  where  two  targets  are  oresent,  but 
with  only  a  small  separation  between  them.  Let  the  first  tarnet  have  an 
intensity  of  unity  and  a  position  at  eauilibrium  of  C(x^iLi),  and  the 
other  target  an  intensity  of  N  and  a  nosition  at  eruilibrium  of 
C(-X2>0),  Then  if  by  annlying  en.  12  to  both  tameto,  the 
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filtered  photocell  sinnal  is 


S(F)  = 


(EjJ;  (x^-Nx^) 
2tt 

0 


■2(N+1)  (£+1; (x^-Nx^) 

TT  ZtT 

0  0 


[17] 


For  fcouilibrium  must  be  zero  (refer  to  eo.  8),  that  is 

A 


2(N+r)pCXj^  -  NX^^  =  0, 


IT 


the  solution  of  which  is 


[18] 



[19] 


This  is  to  be  expected  from  the  souare-law  relationship  between  the 
error  signal  and  the  image  error  (r)  derived  in  Section  2.4*  This 
position  can  be  seen  to  be  stable,  as  eo.  18  shows  that  an  increase  of 
Xj,  with  a  consequent  decrease  in  X2>  will  cause  to  become  oositive, 
whereas  if  x^  becomes  negative  so  does  A  positive  error  signal  will 
cause  the  reticle  to  move  to  the  right  and  hence  decrease  the  error  (xi) 
between  the  reticle  centre  and  a  taraet  image  lying  to  the  right  of 
centre . 
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2 • 6  Lxtended  Target 

First  consider  the  case  where  two  targets,  at  the  same  distance 
from  the  centre,  are  separated  by  the  width  of  one  sector,  as  shown  in 
Fig.  7(a),  which  has  been  illustrated  to  show  the  instant  at  which  the 
reticle  and  spatial  references  are  coincident.  The  waveforms  resulting 
from  the  two  individual  images  are  shown  in  Fig.  7(b),  and  the  total 
signal  in  Fig.  7(c),  where  and  are  the  arguments  of  the  nositions 
of  imaoes  A  and  B. 

An  analysis  of  the  waveform  of  Fig.  7(c),  which  is  given  in 
Appendix  A,  shows  that 

2tt  2ti 

sin  pa  •  fCA+B)  da  =  ^  cos  pa  •  £CA+B)  da  =  0  .  [20] 

and  thus  a^^  ”  ^12  “  seen,  from  eo.  8,  that  because  of  this 

the  error  signal  will  also  be  zero.  It  follows  that  if  the  target  image 
was  in  the  form  of  an  arc,  covering  more  than  one  sector,  the  carrie  " 
signals  from  all  points  separated  by  exactly  one  sector's  width  would 
nullify  one  another.  If  the  arc  was  to  cover  exactly  an  even  number  of 
sectors,  the  carrier  would  be  lost,  and  the  error  signal  would  be  zero. 

The  signals  from  an  arc  image  will  now  be  examined.  An  arc  is 
chosen  as  an  image  shape,  in  preference  say  to  a  straight  line,  onlv 
because  of  its  mathematical  simplicity  in  this  introductory  analysis. 
It  suffers  from  the  disadvantage  that  nothing  can  be  concluded  about  the 
stability  of  equilibrium  positions,  as  once  the  image  moves  it  is  no 
longer  an  arc  with  the  same  centre,  and  hence  the  eauations  employed  in 
the  analysis  become  invalid. 
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FIGURE  7  -  Waveforms  from  two  images 
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In  Fin.  C(a),  the  target  image  is  shown  by  a  heavy  are  from 
P(l,-A)  to  P(1,a).  To  derive  the  error  signals  resulting  f ron  this 
image,  the  filtered  signal  for  an  image  of  unit  intensity  situated  at 
P(1,0)  will  first  be  calculated.  Its  nhotocell  signal  is  shown  in 
Fig.  6(b),  and  hence  (see  also  Appendix  B) 


S(F) 


u 


12 

0 


u 


13 

0 


where 


=  7.596/ll7r, 
“l2  “ 

=  7.596/1371. 


[21] 


If  the  intensity  of  the  arc  is  uniform,  and  its  total  radiance  is  unity, 
then  for  a  segment  at  P(1,X)  of  length  6\ 



2A 


l^cos 

nx 

U12COS 

12X 

“13 

IIX 

12X 

'‘is 

[22] 


In  the  limit 


S(F)  *  6X-KiX-^0, 


-A 

A 


U 


11  Sin  llA 
U 


12  sin  12A 
12 


u 


13  sin  13A 
13 

0 


[23] 


FIGURE  8  -  Extended  target 
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Ihe  carrier  amplitude  of  Llie  arc  is 

a^2  =  Csin  12A)/127tA,  [24] 

and  the  sum  of  the  side-band  a:iDiitudes  ir 

7,596  /sin  llA  ,  sin  13a\ 

s  -  —  (-7^  *  ^7-; . 

These  are  plotted  in  Fin.  9.  As  the  peaks  of  both  the  carrier  and  the 
sum  of  the  side-bands  decrease  raoidlv  with  arc  length,  their 
contribution  to  the  total  error  signal,  in  the  presence  of  a  point 
source  of  similar  intensity,  will  be  very  small. 

Further,  an  arc  of  total  length  sav  tt  would  be  in  enuilihrium,  as 
the  carrier  amplitude  is  then  zero,  A  point  image  at  the  centre  of  the 
reticle  is  also  in  equilibrium,  and  eo.  12  shows  that  this  is  because 
the  side-band  amplitudes  are  zero.  However,  if  both  the  arc  and  point 
images  are  present  at  the  same  time,  the  seeker  is  no  longer  in 
equilibrium,  as  the  side-bands  from  the  arc  and  the  carrier  from  the 
point  will  generate  an  error  signal  and  cause  a  displacement  of  the 
seeker . 


I 
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3.0  RESULTS 


3.1  General 


The  preceding  analyses  have  served  either  to  illustrate  the 
functionina  of  the  seeker  or  to  provide  a  means  of  validating  the 
computer  orogramme,  and  as  such  have  investigated  special  cases  or  used 
approximations.  The  rest  of  the  study,  which  will  be  general,  will 
reauire  the  use  of  a  diaital  computer.  The  simulation  model  was 
programmed  in  APL,  and  run  on  the  Xerox  560  digital  computer.  The 
programmes  are  listed  in  Appendix  C. 

3.2  Error  Curves 


The  error  curve  (Fig,  10(a))  gives  the  error  signal  as  a  function 
of  the  distance  of  a  point  imaoe  from  the  centre  of  the  reticle.  It  is 
to  be  noted  that  the  curve  gives  an  overall  aopearance  of  linearitv,  as 
was  suggested  by  the  reasoning  in  Section  2.2.  The  number  of  neaks  on 
the  curve  corresponds  to  the  number  of  changes  of  snoke  interceptions 
with  the  phasing  sector  as  r  varies  from  0  to  1. 

Figure  10(b),  which  is  plotted  from  the  data  of  Table  I,  shows 
the  portion  of  the  error  curve  for  small  values  of  r.  In  Section  2. A 
the  error  signal  was  derived  as  being  equal  to  6r^  for  values  of  r  less 
than  0.08;  this  verifies  that  the  oart  of  the  error  curve  for  small 
values  of  r  is  correct. 


4 


I 
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5 . 3  Equilibrium  Position  for  Two  Point  Targets 

In  this  test,  two  point  target  images  are  placed  at  a  given 
distance  apart  and,  after  eauilibrium  has  been  reached,  their  distances 
from  the  reticle  centre  are  noted.  The  two  images  are  situated  at 
CCx^jO)  and  C(x2»0),  with  intensities  of  1  and  N.  The  images  are  a 
fixed  distance  apart  (d),  where  d  =  Xj^  -  X2. 

From  Fig.  11  and  Table  II  it  can  be  seen  for  d  =  O.DOl  that 

log  N  Ix^/Xjl  =  i  log  N, 

or  1  [ 26 ] 

|xj/x2l  =  N  . 

This  result  is  in  agreement  with  eo.  19. 

For  larger  d  the  value  of  N|x^/xj^|  in  rinse  to  unitv,  that  is 
if  neither  image  is  close  to  the  centre  of  the  reticle,  so 
that  a  weighted  mean  can  be  applied  in  this  case  to  calculate  the 
equilibrium  position.  Once  either  of  the  imaaes  is  close  to  the  reticle 
centre,  the  signals  from  this  image  will  no  longer  be  linearly 
proportional  to  its  distance  from  the  centre;  nevertheless,  the  slopes 
of  the  curves  for  large  N  remains  one  half. 
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N 

FIGURE  11  -  Two  point  targets 
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TABLC  II 


oint  imanes 


» 

X 

2 

lx  /X  ] 

1  j/  2' 

d  =  0.001 

i.no  .oonso 

-  .00050 

1 . 00 

1 . 00 

1.5S  .0005G 

-  ,00044 

1  .26 

1 . 26 

2,51  .00061 

-  .00039 

1 ,58 

1.59 

3. 93  .00067 

-  .00033 

1 . 99 

2,00 

6.31  . 00072 

-  ,00028 

2.51 

2.  51 

10.00  . 00076 

-  .00024 

3.  16 

3.15 

1 5 . 8S  , 00030 

-  .00020 

3 .98 

3.98 

25. 12  .00033 

-  .00017 

5.01 

5.01 

39.81  . 00036 

-  .00014 

6,  31 

6,31 

63.10  .00039 

-  .0001  1 

7 . 94 

7.94 

100.00  . 00091 

-  .00009 

10.00 

10.00 

d  =  0.4 

*1.00  . 200 

-  .200 

1 . 00 

1 .00 

1.26  .225 

-.175 

1 .28 

.98 

*1.53  .249 

-  .  151 

1  .  64 

.  97 

2.00  .269 

-  ,  131 

2.07 

.  97 

2.51  .285 

-.115 

2.48 

1 .01 

3. 16  .293 

-  .  102 

2 . 92 

1 .08 

■'  3.93  ,309 

-  .091 

1 .  on 

3 . 93 

5,01  , 3  1  S 

-  .08? 

3.8  8 

1  .  29 

6.31  ,325 

-  .074 

4 . 44 

1.42 

10.00  , 34C 

-  .060 

5 .72 

1.75 

15.85  ,352 

-  .04  8 

7.3  3 

2.15 

25.12  .361 

-.039 

9.31 

2.70 

39.81  ,369 

-.031 

11.78 

3.38 

63.10  . 375 

-  .025 

15.00 

4.21 

100.00  .330 

-  .020 

18.80 

5.  32 

d  =  1.0 

1.00  .500 

-  .500 

1  .00 

1.00 

1.26  . 576 

-  .424 

1  .  36 

.  93 

1 . 53  .613 

-.387 

1  .  59 

1  .  00 

2.00  .649 

-  .  351 

1.85 

1  .03 

2.51  .693 

-  .307 

2. 26 

1.11 

3,16  .781 

-.219 

3 . 56 

.  89 

3.93  .817 

-  .  1R3 

4,45 

.  89 

5.01  .843 

-  .  157 

5.  37 

.93 

6.31  .865 

-.135 

6,41 

.  98 

10.00  .897 

-  .  103 

8 , 67 

1.15 

15.85  .919 

-  -081 

11.29 

1.40 

25.12  .935 

-  .065 

14.48 

1,73 

39.81  .949 

-  .  051 

18.49 

2.15 

63.10  . 959 

-.041 

23.57 

2.68 

100.00  .  96  8 

-.032 

30.06 

3.33 
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3,4  Representation  of  Line  Targets 


For  various  numbers  of  points  formina  the  line-taroet  imaqe,  the 
values  of  ^23  calculated  and  oresented  in  Table  III  and 
Fig  .  12,  where  the  line  is  0.8  long  and  positioned  symmetrically  about 
the  X  axis  at  x  -  0.5.  It  was  concluded  from  the  figure  that  little 
error  should  occur  by  employing  21  points. 


To  verify  that  this  number  of  points  is  sufficient  to  represent  a 
line  imaqe  anywhere  on  the  X  axis,  the  amplitudes  of  the  carrier  and  the 
sum  of  the  side-band  amplitudes  were  plotted  for  the  same  line  imane, 
and  its  position  along  the  X  axis  was  varied  from  -1.0  to  1.0.  The 
graph  is  given  in  Fig.  13  for  a  line  image  represented  by  21  noints  and 
in  Fig.  14  for  41  points.  The  only  discernable  difference  between  the 
two  curves  is  in  the  peaks  of  the  carrier  immediately  on  either  side  of 
the  ordinate. 


TABLE  III 

Number  of  points  to  form  a  line  image 

Image  length  =0.8  (Symmetrical  about  X  axis), 
Position  of  image;  x  =  0.5  . 


MBER  OF 
INTS 

*12 

*13 

1 

.1099 

-  .5305 

.0929 

3 

.0816 

- .4953 

.1065 

7 

.0232 

-.  0882 

.0113 

11 

.0219 

-.0821 

.0100 

21 

.0213 

-.0/92 

.0094 

31 

.0212 

-.0786 

.  0093 

51 

.0211 

-.0/83 

.0093 

101 

.0211 

- . U7H  2 

.  0  0  9  2 


NUMBER  OF  POINTS 


FIGURE  12  -  Number  of  ooints  to  represent  a  line  imarje 
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3. 5  Cnuilibrium  Position  in  the  Presence  of  a  Line  and  a  Point  Target 

The  difference  between  the  signals  from  the  representation  of  a 
straight-line  image  by  21  and  41  points,  as  shown  in  Figs.  13  and  14,  is 
small,  and  hence  the  following  computations  have  been  done  with  21 
points;  the  effect  of  the  difference,  if  any,  will  be  discussed  at  the 
end  of  this  section. 


In  Figs.  15(a)-17(a)  are  shown  the  values  a,o  and  af.(=  a-,  + 


12 


11 


for  a  point  imaqe  of  intensity  N  and  a  line  imane  of  unit  intensity 
and  of  length  0.8  centred  about  the  X  axis.  It  is  now  required  to 
derive  the  points  of  equilibrium  from  these  curves  for  a  qiven 
separation  of  the  two  imaqes.  The  separation  will  be  taken  to  be  0.5, 
with  the  point  image  lying  to  the  left  nf  the  line  imane.  For 
equilibrium  to  occur  the  error  sinnal  must  be  zero,  that  is 



[27] 


where 


Xp  is  the  position  of  the  point  imaqe, 

X|^  is  the  position  of  the  line  inane, 

a.(x.  )  is  the  value  of  a  .  at  the  point  x. 
J  K  J  K 


on  the  k  curve, 


and 
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FIGURE  15  -  Line  and  point  images; 

separation  =  0.5,  N  =  1 
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L ine  and  point  images ; 
separation  =  0.5,  N  =  0.3 
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,015 


-  Line  and  point  images; 
separation  =  0.5,  N  =  0 
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Thus  the  equilibrium  position  can  be  given  by  the  solutions  of  both 



ag(Xj^  -  0.5)  +  ag(Xj^)  =  0  . 


[28] 


These  equations  could  be  solved  bv  nlottinq  the  error  signal,  and 
noting  the  zero  crossings.  A  positive  slope  at  these  oositions  would 
show  stability,  as  a  movement  of  the  reticle  centre  to  the  left  would  be 
accompanied  by  a  positive  error  signal,  which  would  tend  to  return  the 
reticle  to  its  equilibrium  position.  Although  the  error  signal  can  he 
very  easily  calculated,  the  value  of  RMS  is  involved,  for  which  a^l’  ^]2 
and  a^^  required.  A  more  direct  approach  is  to  use  the  data  from 
the  curves  of  Figs.  15(a)-17(a),  where  only  the  ai-,'s  and  a,-'s  are 


12 


given.  Let 


Ev  =  e'/RMS^  , 


[29] 


2  ' 

Then  because  RMS  cannot  be  negative,  must  have  the  same  zero 


crossings  as  and  the  signs  of  the  slones  at  these  nositions  must  be 

the  same  for  Ey  and  f'.  Thus,  the  eouilibrium  nositions  are  determined 

A  X 

I 

by  plotting  where 



[30] 


In  Figs.  15(b),  16(b)  and  17(b)  this  curve  is  shown  for  N  =  l.n,  0.5  and 
0.1  respectively. 

It  can  be  seen  for  N  =  1.0  that  there  is  one  eauilihrium 
position,  that  it  is  stable,  and  occurs  for  Xj^  =  0.^5,  Xp  =  -0.07.  For 
N  =  0.3,  there  are  three  ec/ui  1  ihr ium  nositions,  two  of  which  are  stahlp. 
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For  N  =  0.1,  there  are  15  eauilibriun  oositions  of  which  8  are  stable. 
The  figures  show  that  if  a  41-noint  representation  had  been  emoloved 
instead  of  a  21-point  one,  the  larger  carrier  peak  to  the  right  of  the 
ordinate  might  have  caused  two  more  equilibrium  positions  to  have 
occurred  at  x^s  between  0  and  0.1,  for  N  =  0.3. 
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A . G  DISCUSSION 

In  practice,  a  small  nortion  of  the  centre  of  the  reticle  may  be 
blacked  out  to  prevent  oscillation  of  the  seeker  when  the  error  is 
small.  This  would  change  the  shape  of  the  two  point  target  equilibrium 
curves  in  the  region  where  one  of  the  imaaes  enters  this  zero 
transmission  zone. 

Multiple  statically  stable  equilibrium  positions  were  found  under 
certain  conditions  when  both  a  point  and  a  line  image  were  present  on 
the  reticle.  Under  dynamic  conditions  some  of  these  positions  may  not 
be  stable. 


The  filters  employed  extracted  the  carrier  and  its  immediate 
side-bands  with  no  attenuation  or  phase  shift.  The  simulation  model 
could  be  readily  adapted  to  filters  requiring  attenuation  and  phase 
shift  at  particular  harmonics. 

The  only  extended  target  examined  was  a  line.  However,  areal 
targets,  such  as  illustrated  in  Fig.  18,  could  be  studied,  but  because 
of  the  large  number  of  points  required  to  define  an  areal  target  and  the 
conseouent  large  number  of  computer  operations  entailed,  the  programming 
would  probably  have  to  be  performed  in  a  language  other  than  APL. 

The  reticle  employed  is  the  *sunburst*  reticle  with  straight 
spokes.  Other  reticles  could  be  employed  with  very  small  changes  to  the 
simulation.  One  such  reticle  is  illustrated  in  Fig.  19,  in  which  the 
spiral  spokes  improve  the  spatial  filtering. 


FIGURE  18  -  Distributed  target 


FIGURE  19  -  Spiralled  reticle 


I  i 
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5.0  CONCLUSION 

A  simulation  model  of  a  generic  seeker  employing  an  AM  reticle 
was  developed,  from  which  a  digital  programme  was  prepared  to  calculate 
the  seeker  response  to  various  stimuli.  It  was  demonstrated  that  the 
eguilibrium  position  of  two  point  images  can  be  described  by  a  weighted 
mean  only  if  the  separation  and  intensity  ratio  is  such  that  neither 
image,  at  equilibrium,  is  nearer  the  reticle  centre  than  0.1  of  the 
reticle  radius.  This  occurs  because  the  system  is  linear  except  for 
images  near  the  centre  of  the  reticle. 

It  was  shewn  that,  because  of  phase  cancellations,  an  extended 
image  contributes  less  to  the  error  signal  than  a  point  of  the  same 
total  intensity;  hence  a  weiahted  mean  cannot  be  used  for  determininn 
the  equilibrium  position  of  an  extended  and  a  ooint  imaoe.  The 
programme  disclosed  that  with  both  a  line  and  a  point  imaqe  present  on 
the  reticle,  for  certain  ranges  of  intensity  ratio,  there  could  he  many 
static  equilibrium  positions. 

The  model  could  be  modified  to  accommodate  the  given 
characteristics  of  a  particular  AM  reticle  seeker,  such  as  the  geometry 
of  the  reticle,  the  amplitude  and  phase  responses  of  the  carrier  filter, 
and  the  design  of  the  automatic  gain  control.  Areal  targets  of 
specified  shapes  and  gradations  in  intensity  could  also  be  simulated. 


I 
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APPENDIX  A 


Derivation  of  Equation  20 


From  Fig.  7 


f(B,a)  =  f(A,  a+7r/12). 


[Al] 


where 


a  =  u)  t  . 
s 


It  follows  that 


2tt  2Tr 

J  f(B,o)sin  12o  do  =  /* f (A,o+ii/12)sin  12o  do, 

/V 


Theref  ore 


2tt 


-  J  f  (A,a)sin(12a  -  tt)  da, 
0 


2tt 


f(A,a)sin  12a  da. 


[A2] 


2it  2it 

y  f(A-»-B)sin  12o  do  =  y  ff (A)  +  f(B)]sin  12a  da 


0 


»  0, 


[A3] 


Similarly 


2ir 

y  f (A+B) 
0 


cos  12a  da 


[A4] 


( 
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APPENDIX  B 


Derivation  of  the  Constant  in  Equation  21 


From  Fiq.  8 



2t\ 


(photocell  signal)sin  pa  da. 


1  ^ 

-  -r.  ^  (-1)  [cos'rr^pTT  -  cos 


i=0 


(Zli) 

12 



[Bl] 


TT, 


=  —  (1-COS  fr)  ^  cos(7+2j)^ 
np 


"  ^  (l+2cos|^). 


Therefore 


=  7.596/llTr, 
ai2  *  -  Vir, 


a 


13 


[B2] 


=  7.596/13Tr. 


r 
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APPENDIX  C 


Listing  of  Computer  Programmes 


Function  QL2 


VOUT*-QL2\I-,CiA\Pi  Q;  5;  T;  A’VS 

[1]  I-^l 

[2]  P-'-C’-^O 

[3]  (Ar;ircn>‘A’A:[i]  )  +  yyr.r]xyy[i]  )*. s 

[4]  -*-(2'>l  )/A'r 

[53  Axcrj.yyci] 

[6]  A-^l  2®.o/7CxFl 

[7]  P®-P+(C[  l;3x  A[  2;  ])-C[  2  ;3xA[l;] 

[  s]  l;3xA[l;])  +  PC2;]xA[2;1 

[G]  .V2’:-(I  =  pyA')/py 

[10]  Z-«-I+l 

[11] 

[12]  P-V:P’>f5®-(+/5x5®-P.i)*.  5 

[13]  P-PtR:^S 

[14] 

[15]  .9®-+/"l+Pxl(|>P 

[16]  5®-5+  +  /~l+Cxl(t>a 

[17]  3-*-*-/“l+tfxl4»P 

[18]  r«-3'-'*-/“l*Pxl4>i? 

[19]  OiJT^l-S},T 
V 


Use:  Calculates  E^  and  Ey  for  any  number  of  point  images  defined  by  the 
vector  XX,  YY  and  ZZ,  where  XX(I),  YY(I)  and  ZZ(I)  correspond  to 
the  position  and  intensity  of  the  i'th  image. 


Calls:  QL,  QLl 
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Function  QLl 


VG-^QLl  RRiYiPiQiCOiEiTERMiN 
C 1  ]  Y*QL  RR 

[2] 

[3] 

[4]  CO-^loPLl2>tCG*(  .  5«iV^24)*(  -1+L  .  S>iTERM)*\TERM-^3 

[5]  y-(o,fl[i]=o)*y 

[6] 

C7]  ^♦♦[ll^Cd.SxpP)  .2)pP 

[8]  (^-♦/C2]-/C2]2  \».OB».*GG 

C9]  Gl2il-^Gi2il*CO 
CIO]  <;Cli]-— C[l;] 

[11]  (P>-G*oi2,TERIil)fiGG 
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Use;  Determines  b  (i  =  11,  12,  13)  for  a  given  point  image 


Calls:  QLl 


Called  by;  QL2 


UNCLASSIFIFID 

4A 


Function  QL 


V'l*-QL  RR\k',X\  7.\P\C\  iWiH 
[  1  ]  n-*-  (  +  /  A  a  *  /v  A  )  *  .  5 

[2]  fi-'-A'AL  l]Ai'A(V  /VA'[2] 

[3]  A-'-Cot.Sx/iOxO,  i“l+//-^24 
:4]  A^(l.pA)pA 

L5D  A-hA,[1](  .pA)pO  1 

[6]  A'[1;>A[1;]+(o2)xP-^A[1;]<0 

[7]  A-h(+/P)4)A 

[8]  A[  1 ;  ]-AC  1 ;  ]- (  02)x  P-^-AC  1  ;  l>o2 

[9]  a<^(-  +  /p)4;a 

[10]  y<-(  0,  ( ( i+py)  pi )  .o)\y 

[11]  y[l;l]'^0 

[12]  y[l  ;l  +  py]-^o2 

[13]  r[2;l]-l-y[2;2] 

[14]  y[2  ;i*pr]-^i-y[2  ;"i  +  i  +  py  ] 

[15]  2-^(0. 5)-"2oP 

[16]  P-y[l;] 

[17]  ii-P>Z 

[18] 

[19]  -►(fc'^'=l  )//IO 

[20]  fCt^-l]<^l 

[21]  y«-t/y 

[22]  y[l;l]-^Z 

[23]  AD;Z<-(02)-Z 

[24]  P>-y[l;] 

[25]  i-P<Z 

[26]  CQ-C\Q 

[27]  -(4;i)>pC)/0 

[28]  4'[55]-»-l 

[29]  y-i/y 

[30]  y[l;"ltpy]>H2; 

V 


Use;  Determines  zero~cro8sinp  positions  of  the  photocell  waveform  for 
a  given  uoint  image . 


Called  by:  QLl 
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UNCLASSIFIED 

45 


Function  MP 


VI^.VP  NiAiJiMOliPOi. 

[1]  M01^LINBLi2  3] 

[2]  P01-^PI«[;2  3] 

[3] 

C  4  ]  A-^xO 

[5]  i1S:J-^I+50 

[6]  .(MOICI;  1  ]  +  POl[  1  ]*ff )  xMOl  [ I ;  2  ]+ POl  [ ;  2  ]k /y 

[7]  ■*-iJs200)/ AD 
C8]  I-Kl+l 

[9]  -►/IS 

[10]  ADiI^~A 

V 


Use: 


Calculates  the  data  for  Figs  15(b),  16(b)  and  17(b)  from 
matrices  LINE  and  PIN,  where  LINE  contains  the  values  of  a, 2 
and  ag(X|^)  for  -Isx.sl,  and  PIN  contains  the  values  of  3^2 
and  ac(>‘n)  for  si. 

W  P  P 


the 

(x^) 




... 




UNCLASSIFIED 

A6 


Function  ATAN 


VZ*-X  ATAN  T 
Cl]  -*•(  lff"6<U)/I/l 
12}  X*0 

[3]  UliZ-^~ZOY*X*lE  7 

[4]  *(x<o)/yp 

C5]  ^AB 

[6]  YPiZ-rZ*0\ 

[7]  /flfl;  Z>*-(02)  IZ 

[8]  -►(ZSOl)/0 

[9]  Z<-Z-02 
V 


Use;  Calculates  the  arctangent  and  assigns  to  it  the  appropriate  value 
dependent  on  its  guadrant. 
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